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There  are  also  disadvantages:  (1)  Tests  will  not  reveal  fixed-base  natural  frequencies  if  the  3yBtem  con¬ 
tains  extensive  viscous-type  damping.  This  is  not  a  serious  limitation,  however,  since  in  the  presence  of 
large  amounts  of  energy  dissipation,  natural  freuqendes  do  not  in  any  case  represent  the  significant  dynamic 
parameters.  (2)  Precise  setting  and  alignment  of  the  shaker  is  essential  to  avoid  introducing  unwanted  com¬ 
plex  responses. 

In  summary,  this  paper  presents  a  practical  method  to  determine  the  fixed-base  natural  frequencies 
of  a  subsystem.  If  such  frequencies  exist,  they  are  significant  in  current  methods  of  structure  analysis  and 
design  criteria  involving  dynamic  loading. 
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13  ABSTRACT 


The  fixed-base  natural  frequencies  of  a  system  are  those  natural  frequencies  the  system  would  have  if 
it  were  mounted  on  a  base  of  infinite  mass  and  stiffness.  Hie  main  purpose  of  this  investigation  is  to  find  a 
generalized  analytical  method  for  determining  the  fixed-base  natural  frequencies  of  an  in  situ  or  laboratory 
mechanical  subsystem.  The  analysis  will  be  based  on  an  arbitrary  n -degree-of -freedom  undamped  linear 
time-invariant  system.  The  application  of  this  analytical  method  is  extended  to  cases  with  small  viscous 
damping.  Only  one  shaker  is  needed  to  obtain  the  required  response  measurements  necessary  to  calculate 
the  fixed-base  natural  frequencies  of  the  subsystem.  Gages  to  obtain  these  measurements  should  be  placed 
at  all  supporting  points  of  the  subsystem  being  tested  and  one  additional  point  located  any  other  place  on 
that  subsystem.  The  analytical  work  reveals  that  the  following  advantages  of  this  method  are  realized. 

(1)  No  assumption  other  than  linear  time-invariant  of  the  total  system  has  to  be  imposed  in  the  derivation. 

(2)  No  detailed  physical  properties  of  the  whole  system  have  to  be  known  although,  by  its  very  nature,  the 
mobility  of  any  element  within  the  system  depends  upon  the  physical  properties  of  the  whole  system.  In 
fact,  acquiring  detailed  physical  properties  has  always  been  the  difficulty  in  calculating  the  fixed-base 
natural  frequencies  of  a  real  complex  system.  (3)  It  is  equally  convenient  to  carry  out  the  test  in  the  field 
as  well  as  in  the  laboratory,  provided  that  the  shaker  is  capable  of  generating  enough  mechanical  energy 

to  excite  tne  subsystem.  (4)  When  conditions  allow,  many  subsystems  can  be  measured  simultaneously. 
Consequently,  based  on  the  liberal  choice  of  force  application  points,  savings  in  manpower,  time,  equipment, 
and  money  may  be  effected, 

(Continued  next  page) 
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ABSTRACT 


The  fixed-base  natural  frequencies  of  a  system  are  those  natural  frequencies  the  system 
would  have  if  it  were  mounted  on  a  base  of  infinite  mass  and  stiffness.  The  main  purpose  of 
this  investigation  is  to  find  a  generalized  analytical  method  for  determining  the  fixed-base 
natural  frequencies  of  an  in  situ  or  laboratory  mechanical  subsystem.  The  analysis  will  be 
based  on  an  arbitrary  n-degree-of-freedom  undamped  linear  time-invariant  system.  The  appli¬ 
cation  of  this  analytical  method  is  extended  to  cases  with  small  viscous  damping.  Only  one 
shaker  is  needed  to  obtain  the  required  response  measurements  necessary  to  calculate  the 
fixed-base  natural  frequencies  of  the  subsystem.  Gages  to  obtain  these  measurements  should 
be  placed  at  all  supporting  points  of  the  subsystem  being  tested  and  one  additional  point 
located  any  other  place  on  that  subsystem.  The  analytical  work  reveals  that  the  following 
advantages  of  this  method  are  realized: 

1.  No  assumption  other  than  linear  time-invariant  of  the  total  system  has  to  be 
imposed  in  the  derivation. 

2.  No  detailed  physical  properties  of  the  whole  system  have  to  be  known  although, 
by  its  very  nature,  the  mobility  of  any  element  within  the  system  depends  upon  the  physical 
properties  of  the  whole  system.  In  fact,  acquiring  detailed  physical  properties  has  always 
been  the  difficulty  in  calculating  the  fixed-base  natural  frequencies  of  a  real  complex  system. 

3.  It  is  equally  convenient  to  cany  out  the  test  in  the  field  as  well  as  in  the  laboratory, 
provided  that  the  shaker  is  capable  of  generating  enough  mechanical  energy  to  excite  the 
subsystem. 

4.  When  conditions  allow,  many  subsystems  can  be  measured  simultaneously.  Con¬ 
sequently,  based  on  the  liheral  choice  of  force  application  points,  savings  in  manpower,  time, 
equipment,  and  money  may  be  effected. 

There  are  also  disadvantages: 

1.  Tests  will  not  reveal  fixed-base  natural  frequencies  if  the  system  contains  exten¬ 
sive  viscous-type  damping.  This  is  not  a  serious  limitation,  however,  since  in  the  presence  of 
large  amounts  of  energy  dissipation,  natural  frequencies  do  not  in  any  case  represent  the 
significant  dynamic  parameters. 

2.  Precise  setting  and  alignment  of  the  shaker  is  essential  to  avoid  introducing  un¬ 
wanted  complex  responses. 
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In  summary,  this  paper  presents  a  practical  method  to  determine  the  fixed-base  natural 
frequencies  of  a  subsystem.  If  such  frequencies  exist,  they  are  significant  in  current  methods 
of  structure  analysis  and  design  criteria  involving  dynamic  loading. 
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SYMBOLS 


The  syrr  jls  used  in  this  report  are  defined  as  they  appear  within  the  text.  The  most 
important  ones  are  listed  below  for  reference. 

A  dot  over  a  variable  in  the  text  represents  the  differentiation  of  that  var  able  with  re¬ 
spect  to  time. 

A  superscript  T  or  -1  at  the  right  corner  of  a  matrix  indicates  the  transpose  or  inverse 
of  the  matrix  respectively. 

row  matrix 

{  }  column  matrix 

[  ]  square  matrix 

|  ]  non-square  matrix 

[/]  identity  matrix 

|  I  determinant 

||  ||  absolute  value 

Oy  matrix  element  of  [A  ]  at  ith  row  and  ;th  column 

L{q,q;t)  Lagrangian 

T{q;t)  kinetic  energy 

V(q;t)  potential 

qk  feth  generalized  coordinate 

t  time 

[Af]  mass  matrix 

[  R  ]  damping  matrix 

[JC]  stiffness  matrix 
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frequency  of  excitation 

jth  fixed-base  natural  frequency  of  structure  or  substructure  i 
force  at  location  fe 

amplitude  of  the  sinusoidal  force  Ffe(f) 

Spatial  force  column  matrix  (or  vector) 

relative  phase  angle  of  feth  response  with  respect  to  the  applied  force  Fj(t) 
at  location  j 

amplitude  of  the  feth  displacement  response  due  to  the  applied  force  Fj{t) 
at  location  j 

feth  spatial  displacement  component 
Spatial  displacement  column  matrix  (or  vector) 
flexibility  matrix 

flexibility  matrix  element  ai  ith  row  and  ;th  column 
mobility  matrix 

mobility  matrix  element  at  ith  row  and  ;th  column 
cofactor  of  the  flexibility  matrix  element  Zji 
flexibility  matrix  of  the  equipment  subsystem 
flexibility  matrix  of  the  support  subsystem 
flexibility  matrix  of  the  base  subsystem 
coupling  matrix  between  subsystem  fe  and  subsystem  8 


Resonance  function 


DETERMINATION  OF  FIXED-BASE  NATURAL  FREQUENCY  OF 
MULTIPLE-FOUNDATION  MECHANICAL  SYSTEMS  BY  SHAKE  TEST 


INTRODUCTION 

The  theory  explaining  the  small  oscillation  of  a  simple  mechanical  system  has  been  ex¬ 
perimentally  validated  for  some  time.  The  application  of  the  theory  to  a  real  complex  sys¬ 
tem,  however,  is  limited,  where  precise  information  concerning  the  coupling  mechanisms 
among  the  elements  of  the  system  is  lacking.  Even  when  the  exact  physical  constants  are 
given,  the  numerical  calculations  are  very  cumbersome.  As  a  general  practice,  certain  impor¬ 
tant  constants  can  be  extracted  quantitatively  from  the  experimental  results.  The  shake  test 
is  often  used  to  measure  the  vibrational  responses.  It  was  mentioned  in  the  report  by  Petak 
and  O’Hara  (1)  that  knowledge  of  the  fixed-base  natural  frequencies  of  a  system  is  essential 
in  current  methods  of  treating  dynamic  problems.  Since  Ref.  1,  the  method  of  determina¬ 
tion  of  fixed-base  natural  frequencies  by  shake  test  has  been  investigated  at  NRL.‘  Simple 
cases  were  analyzed  *»y  £etak  Kapian.^2)  for  a  one-dimensional  linear  chain,  and  by-Petak 
and  O’Hara  (1)  for  a  dual-foundation  shipboard  equipment.  Experimental  work  on  a  single 
support  beam  was  reported  by  Remmers  (3).  This  analytical  method  is  proved  to  be  a  gen¬ 
eralization  of  the  aforementioned  special  cases. 

This  report  includes  the  analytical  presentation  of  this  method,  the  physical  interpreta¬ 
tion  of  the  entities  derived,  the  technical  method  suggested  in  taking  and  analyzing  test  data 
and  drawing  conclusions,  and  the  extension  of  its  application  to  damped  systems,  which  is 
of  importance  in  engineering  practice.  Various  computer  simulations  of  a  problem  designed 
to  illustrate  the  developed  method  are  provided.  ’Experimental  confirmation  of  this  method 

on  a  three-support  beam  mounted  on  a  truss-like  frame  will  be  presented  in  a  separate  report.* 

\ 

\ 

The  experimental  techniques  in  principle  are  essentially  the  same  as  described  in  Ref.  2, 
except  that  phase  differences  (or  the  relative  phase)  between  the  applied  force  and  corres¬ 
ponding  responses  at  points  of  measurement  are  necessary  for  a  subsystem  having  three  or 
more  supporting  points.  By  use  of  the  resonance  condition,  the  fixed-base  natural  frequencies 
can  be  calculated  from  the  experimental  data.  In  order  to  eliminate  the  extraneous  frequen 
cies  inherent  in  this  method  of  measurement,  one  additional  shaking  point  on  the  equipment 
subsystem  itself  is  required.  Contrary  to  the  speculation  that  many  shakers  with  variable 
forces  (both  magnitude  and  phase)  ought  to  be  used  in  a  multiple-foundation  system,  this 
method  requires  only  one  shaker  at  a  time.  The  magnitude  of  the  forces  applied  is  irrelevant 
to  the  result,  as  long  as  the  system  is  being  excited  to  a  measurable  level.  However,  the  phase 
differences  between  the  forces  and  the  responses  are  important. 
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THEORETICAL  ANALYSIS  OF  A  LINEAR  TIME-INVARIANT  MECHANICAL  SYSTEM 

It  is  intended  to  develop  a  steady-state  method  for  determining  fixed-base  natural  fre¬ 
quencies  of  an  in  situ  or  laboratory  mechanical  subsystem.  In  this  analytical  treatment,  gen¬ 
eralized  coordinates  q*  in  configuration  space,  Lagrangian  mathematical  formulation,  and 
linear  space  matrix  representation  will  be  used  to  clarify  the  physical  contents  of  the  mathe¬ 
matical  formulation.  Since  it  is  understood  that  the  main  goal  of  this  investigation  is  to  deter¬ 
mine  fixed-base  natural  frequencies,  the  result  is  confined  to  linear  time-invariant  systems, 
or  to  systems  behaving  approximately  so. 


Undamped  System 

Equation  of  Motion— The  Lagrangian  formulation  describing  an  n-degree-of -freedom 
dynamical  system  can  be  found  in  many  books  on  mechanics  (4,  5).  Since  this  formulation 
is  not  our  main  interest  here,  only  the  result  will  be  introduced.  Suppose  L(a,q;t)  denotes 
the  Lagrangian  of  the  system,  T(q;t)  the  kinetic  energy,  and  V(q;t)  the  potential  energy  of 
the  system;  then  by  definition, 


L  =  T(q;t)  -  V(q;t) 

(1) 

T(q;t)  =  |  ^{M]{q} 

(2) 

mt)  =  |  £[*]{<*}. 

(3) 

From  the  variational  principle,  the  equation  of  motion  of  a  free  system  can  be  derived  as 


dai  dL  _ 

dt dqk  dqk 

When  there  are  applied  forces,  the  constitutive  equation  may  be  written  as 


(4) 


d  dL  dL  ... 

dt  dqk  dqk  k 

(5) 

The  explicit  expressions  in  matrix  form  of  Eq.  (4)  and  (5)  are 

[M]{q]  +  [K){q}  -  {0} 

(6) 

[M]{q]  +  [*]{«}  =  {F(t)} 

(7) 

respectively. 

For  the  particular  purpose  of  natural  frequency  determination,  the  applied  forces  F(t) 
on  the  right  side  of  Eq.  (7)  have  to  be  limited  to  being  sinusoidal,  so  as  to  provide  a  single 
input  with  a  well-defined  driving  frequency  u>  to  the  system.  Let 
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F(t)  =  F  smut  .  (8) 

Here  the  phase  of  the  force  F(t)  is  assumed  to  be  zero.  This  assumption  will  not  lose  its  gen¬ 
erality  in  the  analysis,  if  a  relative  phase  is  assigned  to  the  responses.  Then  the  steady-state 
solution  of  the  set  of  n  simultaneous  linear  differential  equations,  Eq.  (7),  has  the  general 
form  of 


{Qkj}  =  {  9*;sin  M  +  <Pkj)}  (9) 

where  <t>kj  =  nn ,  n  =  0, 1, 2,  •  •  •  for  undamped  cases.  Substituting  Eq.  (8)  into  Eq.  (7),  the 
set  of  n  simultaneous  linear  differential  equations  is  reduced  to  a  set  of  n  simultaneous  alge¬ 
braic  equations,  which  has  the  following  general  form 

(-o>2[M]  +  [K]){q}=  {f}  (10) 


where 


-co2liif]  +  [*]  =  [Z]  (11) 

is  defined  as  the  flexibility  matrix  of  the  system,  and  (l/w)[Z)  is  normally  defined  as  me¬ 
chanical  impedance  matrix  (1). 


Linear  Space— In  order  to  simplify  the  mathematical  manipulation  and  help  to  shed 
light  on  the  physical  insight  of  the  formulation,  Eq.  (10)  can  be  written  in  the  following  form: 

[Z]{q}  =  {F}.  (12) 

If  one  views  Eq.  (12)  as  a  mapping  or  transformation  between  two  linear  spaces,  then  [Z]  is 
just  an  operator  which  maps  the  force  space  onto  the  coordinate  space.  Here  [Z]  is  assumed 
to  be  nonsingular;  this  is  true  for  most  real  physical  systems,  and  the  mapping  is  in  a  one-to- 
one  c  jrrespondence.  When  the  inverse  of  [Z]  exists,  one  can  apply  operator  [Z]-1  to  Eq. 
(12)  to  obtain 

[Z]-1[Z]{9}  =  [ZH{f}  ,  (13) 

or 

{<?}=  [zH{f}  , 

where  the  identity  equation 


[Z]-1[Z]  =  [I]  (14) 

is  used.  The  matrix  [Z]-1  is  often  defined  as  the  influence  matrix.  The  mobility  matrix  is 
defined  in  terms  of  the.  influence  matrix  as 
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[*]  ^.[Z]-1.  (15) 

For  the  convenience  of  discussion,  the  linear  spaces  and  {f}  are  divided  into  sub¬ 
spaces.  In  doing  so,  a  few  definitions  are  necessary.  Any  physical  system  may  be  described 
as  three  main  parts:] 

Equipment:  The  subsystem  under  investigation 
Support:]  The  subsystem  supporting  the  equipment 
Base:  The  remaining  part  of  the  total  system  other  than  the  equipment 

and  support. 

Correspondingly,  these  three  parts  can  be  defined  as  three  subspaces  {9*},  {9*}.  {9*},  in  the 
coordinate  space  {<?},  and  three  subspaces  {Fe},  {f*},  {f6}  ,  in  the  force  space{F}.  Then 
Eq.  (12)  can  always  be  arranged  to  have  the  partitioned  form 


"  IZ'1  I  {0} 

w 

{c;}T  |  iz-i  |  {c;}  ■ 

M 

_  {0}  I  {c;rj  (z‘1 

i«6i 

or  by  matrix  multiplication 


[Ze){q*}  -  {Fe}  , 

{C/FM  +  lZs]{qj  +{Cg}{q«’}  ={F*}  , 
{CgfM  +  iZb){qb}  ={F&}. 


(16) 


(17) 

(18) 
(19) 


It  is  apparent  that  there  is  no  direct  coupling  between  the  equipment  and  base  subsystems  be¬ 
cause  of  the  definition  of  the  support  subsystem.  In  the  analysis  as  well  as  tests,  the  three 
subsystems  (or  subspaces)  must  be  defined  without  ambiguity. 


Resonance  Condition— The  resonance  condition  for  the  equipment  subsystem  requires 

[Ze]{qe}  =  0  .  (20) 

Consequently,  this  requires 

{c/}{9s}  =  {f*}  .  (21) 

Theoretically,  there  is  a  definite  solution  for  Eq.  (21)  if  the  rank  of  {C*}equals  the  dimen¬ 
sion  of  { qs }  for  a  given  set  of  {Fe}  ..  But  this  does  not  usually  occur  in  practical  applications. 
However,  the  applied  forces  are  arbitrary  and  under  our  control.  Therefore  one  may  set 


NRL  REPORT  7300 


5 


M  =  W,  that  is,  there  are  no  forces  which  can  be  applied  to  the  equipment,  and  still 
have  the  resonance  condition  satisfied.  Equation  (21)  becomes 


(C*}M  =  0  •  (22) 

Since  {C/}  ^  {o},  the  support-displacement  response  vector  must  vanish  at  the  fixed-base 
natural  frequencies  of  the  equipment  subsystem.  Now  let  us  examine  the  set  of  simultane¬ 
ous  equations,  Eqs.  (17),  (18),  and  (19),  under  the  resonance  condition  discussed  previously. 
The  system  can  be  described  by  the  following  equations: 

[Ze]{q'}  =  0  ,  (23) 

{CffU*}  +  =  {F*}  .  (24) 

[Zb]{qb]  =  {Fb}  .  (25) 

It  can  be  explicitly  shown  that  the  solution  of  Eqs.  (23),  (24),  and  (25)  does  contain  the  in¬ 
formation  we  need  to  know  about  the  fixed-base  natural  frequencies  of  the  equipment  sub¬ 
system.  This  information  is  in  agreement  with  the  computer  simulation  of  a  one-dimensional 
linear  chain  structure  reported  in  Ref.  3.  In  the  present  report,  we  have  shown  analytically 
that  it  is  true  for  all  undamped  linear  time-invariant  systems  in  general,  no  matter  what  the 
forces  are. 


Resonance  Function— A  resonance  function  is  defined  in  this  section  because  of  its  par¬ 
ticular  importance  in  fixed-base  natural-frequency  determination.  It  contains  information 
about  the  resonance  frequencies  of  the  equipment  subsystem.  All  extraneous  natural  fre¬ 
quencies  of  the  base  and  support  subsystems  are  excluded.  The  graphical  representation  of 
the  function  shows  peaks  at  the  fixed-base  natural  frequencies  of  the  equipment  subsystem. 
The  resonance  function  is  also  a  function  of  the  mobility  entities  of  the  mobility  matrix  of 
the  total  system.  These  entities  are  the  characteristic  and  invariant  properties  of  the  total 
physical  system,  and  their  measurability  is  of  key  importance  in  this  investigation. 

Before  we  proceed  to  derive  an  explicit  expression  of  the  resonance  function,  it  seems 
essential  to  recall  the  linear  space  transformation  in  matrix  form.  The  harmonic  solution  has 
shown  that  the  mobility  matrix  is  defined  as  in  Eq.  (15).  The  entities  of  the  mobility  matrix 
can  be  expressed  as 


mkj  =  «  -  (26) 

where  Zjk  is  the  adjoint  of  the  flexibility  matrix  [Z] ,  which  corresponds  to  the  flexibility 
matrix  elemental .  Physically  there  is  a  difference  in  meaning  between  the  mobility  matrix 
and  the  inverse  of  the  flexibility  matrix;  the  mobility  matrix  is  an  operator  to  map  the  gen¬ 
eralized  velocity  vector  onto  the  generalized  force  space,  while  the  inverse  flexibility  matrix 
is  an  operator  to  map  the  generalized  coordinate  vector  onto  the  generalized  force  space.  In 
our  particular  case,  generalized  coordinates  and  generalized  velocities  are  related  by  a  scalar 
multiplier  w,  the  frequency  of  excitation.  It  is  immaterial  which  one  is  used  in  this  analysis. 
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However,  in  order  to  avoid  ever  increasing  confusion  in  terminology,  the  mobilities  are  used 
throughout  the  derivation. 

We  consider  a  system  composed  of  an  8  degree-of-freedom  equipment  subsystem,  an 
m-e  degree-of-freedom  support  subsystem,  and  an  n-m  degree-of-freedom  base  subsystem. 
The  explicit  form  of  the  set  of  simultaneous  equations,  Eq.  (13),  describing  the  total  n 
degree-of-freedom  dynamic  system,  is  the  following,  when  there  is  no  force  acting  on  the 
equipment  subsystem:: 

uqi  =  mx«+iF£+1  +  -  -•  +  mlmFm  +  •••  +  mlnFn  , 
ux?2  =  m2  8+l*8+l  +  ■  ■  ■  +  m2mFm  +  •  •  ■  +  m2  n^n  » 


~  mt  e+i^e+i  +  -  •  •  +  m t  mFm  +  —  +  mg  nFn  , 

<*>98+ 1  =  "*8+1  s+i^c+i  +  •••  +  mc+imFm  +  ...  +  me+i  nFn  , 


u9m  mm  e+l^e+l  +  '  +  mm  m  +  •  •  •  +  tnm  nFn  , 

<*>9m+l  =  mm+l  +  ’  '  *•  +  mm+l  m^m  +  •••  +  mm+i  nFn  , 


<*> >9n  ~  mn  s+iF 8+i  +  •  •  •  +  mn  mFm  +  •  •  •  +  mnnFn  . 

In  applying  the  resonance  condition,  we  single  out  the  part  containing  the  generalized  coor¬ 
dinates  of  the  support  subsystem  from  Eq.  (27)  and  set  them  equal  to  zero;  that  is. 


9e+l  =  9s+2  =  '  •  '  ~  9m  "  0  , 


(28) 


or 


™e+l  e+l^e+i  +  •  •  •  +  me+i  mFm  +  •  -  -  +  me+i  nFn  -  0  , 
^8+2  8+l^e+i  +  +  rnt+2mFm  +  •••  +  me+2  „F„  =  0  , 


(29) 


mm  8+1  ^£+1  +  •  •  •  +  mm  mFm  +  •••••+  mm  nFn  -  0  . 


Equation  (29)  is  a  set  of  simultaneous  homogeneous  algebriac  equations.  In  general  n  -  C  > 
m  -  t ;  that  means  an  infinite  number  of  solutions  can  be  developed  for  various  given  forces. 
In  other  words,  the  solution  is  indefinite.  However,  once  again  one  may  take  advantage  of 
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the  fact  that  the  applied  forces  are  under  our  control  during  the  measurement.  If  one  limits 
the  total  number  of  forces  applied,  one  at  a  time,  to  be  equal  to  the  degree  of  freedom  of  the 
support  subsystem  m  —  £,  then  a  nontrivial  solution  exists.  This  solution  requires  that  the 
determinant,  consisting  of  the  mobility  entities  associated  with  the  applied  forces,  vanishes. 
For  example,  if  all  the  forces  are  applied  at  the  support  points,  then 

I  me+i  t+i  *  •  •  m8+i  m  I 


mg+ 2  o+i  -  •  *  mt+2 


mm  s+l  * '  ‘  mm 


-  n(u)  . 


It  is  equally  true  for  any  other  combination  of  m  -  t  forces.  In  other  words,  the  forces  can 
be  applied  on  the  support  subsystem,  on  the  base  subsystem,  or  on  the  combination  of  both.. 
This  liberal  choice  of  application  points  offers  great  flexibility  in  application  of  this  particular 
method. 

It  is  important  to  show  explicitly  that  Eq.  (30)  does  contain  information  about  the  fixed- 
base  natural  frequencies  of  the  equipment  subsystem.  Recall  the  matrix  Eq.  (14).  Hie  flexi¬ 
bility  matrix  and  its  inverse  an  commutative.  The  matrices  in  Eq.  (14)  may  be  partitioned 
and  written  in  the  following  form: 


{c/f  |  [Z‘l 
{0} 


j»}_ 

{CJ} 


mll/w 


•  •  j  •-  •-  *  I  mnn/w 


'!°!° 
"7“T“ 
o  i  -j  o 

— t-  +  _ 
0  I  0  I  / 


If  one  replaces  the  inverse  flexibility  matrix  by 


[/]  '  mi  l+1/w  •  «•  |  {0} 

_  _  _J_  _  _  _ _ _ _  J _ _ 

W[__]  __  jw 

M  ;  •••  ;[;] 


one  has 


{<?}  {o} 


mi  i+i/u.  •  -  ■  !  {0} 


S°1 


{»}_]  r_U")  |{o}|_{o}1  (31) 

{o}  |  = !  {c;'}7 ;  (/]  {(ft}  ! . 


{c%Y\w\ J 


{o}  j  {o} ;  j 


8 


NI  ANDSKCMP 


Taking  the  determinant  of  both  sides  of  Eq.  (31)  and  applying  the  identical  equation 


IUHB]i  =  lUJMlBlI 


(32) 


results  in 


=  [zel[zfc]wm-t , 


or 


= 


l[Z«ll  l|Z»)l 

(Zll  ' 


(33) 


Equation  (33)  shows  that  =  0  not  only  carries  information  about  the  fixed-base  natural 
frequencies  of  the  equipment  subsystem,  that  [Z*]  =  0,  but  also  the  fixed-base  natural  fre¬ 
quencies  of  the  base  subsystem  as  well  as  the  natural  frequencies  of  the  total  system.  These 
extraneous  frequencies  can  be  eliminated  by  constructing  another  determinant  having  an 
order  of  m  -  C  +  1.  It  is  done  physically  by  adding  another  force  applied  at  one  of  the  equip¬ 
ment  points.  If  we  designate  this  determinant  as  f2’(w),  then 


c  •  *  *  mi  m 


me+l  «+l 


=  n'(co)  . 


mm  8  ...  mm  m 


Following  the  same  procedure  in  deriving  J2(w),  one  has 


n'(w) 


i[z*']i  \[zb}\ 
\[z]\ 


(34) 


In  general,  l[Ze  |l  does  not  vanish  at  the  fixed-base  natural  frequencies  of  the  equipment  sub¬ 
system  unless  the  determinant  |[Ze  ]  1  and  the  determinant  |[Ze]l  do  have  common  factors. 
This  will  be  discussed  later. 


Now  we  are  in  position  to  construct  the  resonance  function  in  explicit  form. 

This  is  defined  as 


v('Hfy.w) 


Sl'(w) 

ilze']i| 

f2(w) 

llZ'Hl 

(35) 


Here,,  the  double  line  designates  absolute  value.  The  use  of  the  absolute  values  in  the  reso¬ 
nance  function  evaluation  essentially  converts  the  zero-crossing  points  of  function  l[Zc]\  =  0 
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to  minima.  For  a  physical  system,  \[Ze]\  is  finite  for  finite  w.  The  singularities  of 
corresponding  to  l[Ze]l  =  0  gives  the  fixed-base  natural  frequencies  of  the  equipment  sub¬ 
system. 

As  for  the  identification  of  the  detected  resonance  frequencies  concerned,  the  extrane¬ 
ous  frequencies  of  the  base  subsystem  and  the  total  system  are  eliminated.  As  pointed  out 
before,  when  there  are  decoupled  subsystems  within  the  equipment  subsystem,  ambiguity 
arises  due  to  the  fact  that  \[Ze  ]\  and  \[Ze)l  have  common  factors.  In  turn,  the  resonance 
function  does  not  reveal  all  the  resonance  frequencies  as  expected.  The  remedy 

for  such  ambiguities  is  to  take  one  or  two  more  force  application  points  on  the  equipment 
subsystem  for  cross-check  purposes. 

Although  the  derivation  of  the  resonance  function  is  based  on  the  assumption  that  the 
m  -  fi  forces  are  applied  to  the  supporting  points,  the  actual  locations  of  the  application 
points  of  the  forces  are  irrelevant  as  long  as  they  are  kept  away  from  the  equipment  subsys¬ 
tem.  The  aforementioned  fact  is  not  obvious,  because  the  extraneous  frequencies  induced  in 
the  derivation  of  the  functions  fl(w)  and  depend  upon  the  force-application  points. 
However,  those  induced  extraneous  frequencies  are  common  factors  of  fi(w)  and  Q'(cj).  As 
a  consequence,  the  resonance  function  remains  unchanged.  A  mathematical  proof  for  this 
liberal  choice  of  force  application  points  will  be  given  in  the  next  paragraph. 

To  prove  the  r*bove  statement  of  liberal  choice  of  force-application  points,  we  refer  to 
Eq.  (31),  in  which  the  matrix  replaced  the  inverse  flexibility  matrix  having  its  physical  sig¬ 
nificance.  In  that  matrix,  one  noticed  that  the  mobility  entities  were  replaced  by  unity  or 
zero,  except  those  lying  in  the  m  -  £  columns.  Those  columns  correspond  to  the  locations 
where  forces  are  applied  and  include  the  eigenvalue  function  ft(w)  for  the  resonance  condi¬ 
tion  of  the  equipment  subsystem  under  investigation.  In  this  particular  case,  the  resonance 
condition  is 


9s+l  =  Qs+2  =  • ' '  Qm  =  0 

with  forces  applied  at  the  supporting  points  qg+i ,  Qt+2  '  Qm-  Now  suppose  one  of  the  ap¬ 
plication  points,  say  qw,  is  to  be  changed  to  a  point  on  the  base  subsystem,  say  at  qv,  where 
C  <  w<  m,  and  m<  v<  n,  the  corresponding  matrix  change  is  as  follows: 


!  ml  E+l /w 

*-  *  '■  Tn  i  w/u 

•  ’  ml  ml w  ! 

Mi 

1 

1  {o} 

1 

— 

1 

jo)j _ 

i 

• 

:  {o} 

i 

(oil 

!  i/] 

i  w,,  c+l/tj 

I 

mn  m/w  ' 

frin  w/lo  1 

1 
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to 


— 

[/] 

ml  «+l/u> 

1 

■mi  w-i/u .  •  *  *  ml  w*ilw  •  •  •  mi  m/u  j  0  •  •  ■  mi  v/lJ  •  -  •  0 

1 

_  __ 

-1- 

{0} 

1 

i _ 

:  1 

4- 

1 

1 

1  / 

!  \ 

{0} 

1 

1 

.1  ‘l 

5  mv  v/u> 

.  mn  8+l/a> 

mn u>-iiw'  ‘ '  m„  io+i fo  *  *  •  mn  m/w  |  I 

1 

1 

1 

i  1 

_ 

1 

_ 

The  determinant  of  this  new  matrix  is  also  an  eigenfunction  S2(w)  corresponding  to  the  same 
resonance  condition  as  stated  before,  but  with  the  difference  that  one  of  the  applied  forces 
is  changed  to  a  new  location  at  qv  on  the  base  subsystem  instead  of  at  qw  on  the  support 
subsystem.  If  one  makes  the  corresponding  change  in  Eq.  (31)  and  designates  the  new  matrix 
as  IV]  ,  one  will  have  the  following  matrix  equation: 


(36) 


mm  = 


mi 


— i - 

1  0 •  •  •  cs+i  v  *  •  • 0 
1 
I 

e\T 


{ces} 


{0} 


{0} 


0  •  •  •  cm+\  v  •  •  •  0 


_ {0}_ _ 

cs+l  e+i  *  ’  *  c»+l  0~l  0  cg+i  y+i 


'm+l  m+1 


■zb  0  zh 

zm+ 1  d-1  u  zm+ 1  v+1 
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Taking  the  determinants  on  both  sides  of  Eq.  (36),  one  has 


l[zll 


mt+ 1  e+l/u>  •  •  ■  me+l  W-lltom*+l  W+llto  •  '  '  Wg+1  rn/u  mi*l  v/to 


mm  C+1/cj  "  "  ’  mm  uf-l/co  mm  ui+l/u  '  '  *  mm  m/to  mm  v/to 


=  \[Ze]\ 


Cm  (/  Cjp  ui+i  •  •  'Cm  Cm  u+i  ■  *  'Cu  n 

b  b  _6  b 

cm+ 1  v  cm+ 1  m+1  *  ’  "zm+ 1  u-1  zm+l  u+1  *  *  *  2m+l  n 


cu-l  u  2u-l  m+1  ‘  *  ‘2u-l  u-1  2u- 1  u+1  '  ‘  "2u-l  n 

b  b  b  b 

cu+l  v  2u+l  m+1  *  *  ‘2u+l  u-1  2u+l  u+1  "  '  *2u+l  n 


„  ,b  _b  b  b 

cn  v  zn  m+1  zn  w-1  zn  u+1  '  zn  n 


or  simply 

S2(o») 

where  [Zb]  represents  the  matrix 


\[Ze)\ 

\[Z)\ 


\izb)\ 


to 


m-s 


(37) 


(38) 


Gw  v  cw  m+i  •  •  •  Cm  cw  y+i  *  •  •  cw  n 

b  b  b  b 

cm+ 1  v  zm+ 1  m+1  '  *  zm+ 1  u-1  2m+ 1  u+1  *  ’  *  2m+ 1  n 


b  b  b  b 

cu- 1  u  2 u-1  m+1  '  '  *  2 u-1  u-1  2 u-1  u+1  *  *  ’zu-l  n 

b  b  b  b 

cu+ 1  u  2 u+1  m+1  '  "  '  2 u+1  u-1  2 u+1  u+1  ‘  ’  ' 2 u+1  n 


cn  v  zn  m+1  *  zn  u-1  zn  u+1  ''  ''  '  zn  n 
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Hie  indices  w  and  v  a re  rather  arbitrary  within  the  limit  of  8  <  w  <  m,  and  m  <  v  <  n.  This 
procedure  can  be  repeated  to  interchange  the  force-application  points  between  the  support 
subsystem  and  the  base  subsystem  as  desired.  When  each  additional  application  point  of  the 
forces  is  changed  to  another  location,  Eqs.  (37)  and  (38)  show  that  the  difference  in  the 
eigenfunction  ft(u>)  is  limited  to  the  induced  extraneous  frequency  function  [Zb]  alone,  but 
not  to  the  absolute  value  of  the  ratio  of  the  determinant  \[Ze]\  and  the  determinant  I [Z ]i. 

This  means  that  such  exchange  of  force-application  points  will  not  alter  the  information  of 
the  fixed-base  natural  frequencies  of  the  equipment  subsystem.  Nevertheless,  the  solutions  of 

\[Zb]\  =  0  and  \[Zb'}\  =  0, 

in  general,  will  not  be  identical.  Therefore  the  remaining  task  is  to  show  that  the  extraneous 
frequency  function,  \[Zb]\  =  0,  will  remain  unchanged  if  an  additional  force  is  applied  to  the 
equipment  subsystem.  Attention  is  called  to  the  fact  that  the  flexibility  matrices  [Ze]  and 
[Zb]  of  the  subsystems  are  merely  artificial  designations  for  the  convenience  of  discussion. 
They  a  e  physically  equivalent,  in  reality;  choice  of  one  or  the  other  is  dependent  on  which 
part  of  he  system  is  of  interest  in  the  particular  investigation.  Because  of  this  equivalence 
betwee  .  the  equipment  and  base  subsystems,  one  may  conclude,  without  repeating  the  mathe¬ 
matical  manipulation,  that  one  additional  shaking  point  on  the  equipment  subsystem  will  not 
change  the  determinant  \[Zb]\,  but  will  change  the  determinant  \[Ze\\.  This  is  the  property 
being  used  in  constructing  the  resonance  function  ^(m^-; w)  to  assure  that  all  the  induced  ex¬ 
traneous  frequencies  are  eliminated. 


Damped  System 

The  computer  simulations  in  the  “Example  Problem”  demonstrate  the  soundness  of  the 
developed  fixed-base  natural-frequency  determination  method  for  a  one-dimensional  multiple- 
foundation  undamped  system.  It  is  of  both  theoretical  and  practical  interest  to  ask  whether 
this  method  can  be  used  in  studying  damped  systems.  In  order  to  answer  this  question,  an 
analytical  approach  is  necessary  for  two  reasons:  first,  the  extension  of  this  method  to  apply 
to  damped  systems  is  not  trivial;  second,  damping  does  exist  in  all  physical  systems.  In  other 
words,  this  analysis  is  intended  to  show  the  applicability  and  limitation  of  the  fixed-base 
natural-frequency  determination  method  on  damped  systems  on  the  one  hand,  and  the  in¬ 
trinsic  properties  of  the  damped  system  on  the  other.  Either  of  the  aforementioned  interests 
requires  proper  understanding  of  the  experimental  results.  Such  understanding  cannot  be 
achieved  merely  from  the  measured  dynamic  responses.  The  exact  manner  of  acquiring  the 
test  data  sometimes  is  the  determining  factor  to  accomplish  the  purpose.  It  is  known  that 
the  shake  test  has  been  used  to  measure  impedance,  mobility,  and  natural  frequency  of  a 
mechanical  system  directly  or  indirectly  through  response  measurements.  Those  quantities 
are  immensely  useful  in  structure  analysis  and  design  criteria  where  dynamic  loadings  are 
involved.  The  purpose  of  this  analysis  is  to  provide  a  better  understanding  of  real  complex 
systems,  so  that  use  of  test  data  will  be  more  confident. 


Equation  of  Motion— A  damped  system  involves  dissipation  of  energy.;  In  numer¬ 
ous  experimentally  justified  cases,  the  dissipative  forces  in  a  mechanical  system  are  indeed 
proportional  to  velocities  and  referred  to  as  viscous  damping.  For  such  a  system,  one  may 
use  the  Rayleigh  dissipation  function,  which  is  defined  as 
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R  =  09) 

where  ry  are  the  elements  of  the  damping  matrix  [12  ],  which  are  real  and  positive. 

The  corresponding  constitutive  equation  of  motion  in  the  Lagrangian  formulation 
becomes 


or,  written  in  explicit  matrix  representation  of  the  governing  differential  equations, 

lM]{q}  +  [J2]{q}  -r  [K]{q}  =  {f}  .  (41) 

A  great  number  of  structures  also  exhibit  energy  loss  during  motion,  with  the  dissipa¬ 
tive  forces  proportional  to  displacements  known  as  coulomb  damping.  It  is  equally  possible 
that  energy  losses  due  to  dissipative  forces  proportional  to  acceleration  may  have  a  sizable  in¬ 
fluence  on  the  system  in  motion.  However,  these  latter  cases  may  be  described  by  Eq.  (41), 
if  we  regard  the  mass  matrix  [M  ]  and  the  stiffness  matrix  [K ]  as  effective  mass  matrix  and 
effective  stiffness  matrix  accordingly.  Therefore,  Eq.  (41)  may  be  used  to  study  damped 
linear  time-invariant  systems  in  general. 


Closed  Form  Mathematical  Solution— In  an  idealized  shake  test,  a  shaker  providing  a  real 
sinusoidal  force  is  assumed  in  a  mathematical  model.  Such  a  force  may  be  expressed  as 

F(t)  -  Fsinwf. 


The  steady-state  solution  of  the  governing  differential  equation,  Eq.  (41),  is  supposed  to 
exist,  in  which  the  whole  system  vibrates  at  the  unique  frequency  co  of  excitation.  The  re¬ 
sponses  can  be  expressed  in  the  following  general  form: 


{Qk}  =  {?fe  sin  {cot  +  <pkj)} 

=  {afe  cos  cot  +  bk  sin  cot }  , 


(42) 


where  ak  ■=  q*  cos  0fe;- ,  b*  =  q*  sin  <j>kj ,  and  in  general  <t>kj¥=  m ,n  =  0, 1,  2, . . .  for  cases  in 
which  viscous  damping  is  involved.  It  is  evident  that  in  a  dpmped  system,  two  components  are 
needed  to  describe  a  single  response,  i.e.,  either  ( qk,<t>kj )  or  >  bk )•  Taking  the  time  deriva¬ 

tives  of  Eq.  (42)  and  substituting  them  into  Eq.  (41),  one  has 

-w2[Af]{a  sin  cut  +  bcoswt}  +  u>[l2]{a  sin  cot  -  b  cos  cot) 

+  [*]{a-  sin  cot  +  b  cos  wt}  =  {f  sin  cot }  . 


(43) 
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Since  sin  wt  and  cos  u>t  are  independent,  Eq.  (43)  may  be  written  as 

[Af]{a}  -  u[R]{b)  +  =  {f}  ,  (44) 

-cj2{M]la}-  a>[i?l{&}  +  [K]{i}  =  {o }  .  (4?> 

By  rearrangement  of  the  terms  in  Eqs.  (44)  and  (45),  and  writing  them  in  a  compact  form, 
we  have 

[X]  -  wTO  |  -wllil 

- 1 - 

w[£]  j  [ K ]  -  u^M) 


Theoretically,  the  problem  is  solved  once  a  set  of  forces  along  with  the  physical  properties 
[M],[X],andlF]  of  the  system  are  given.  A  closed-form  solution  can  be  worked  out  from 
Eq.  (46)  by  any  computer  means.  However,  the  application  of  Eq.  (46)  in  practical  measure¬ 
ment  is  not  trivial.  With  this  closed-form  solution  in  mind,  an  approximate  solution  for  a 
shake  test  will  be  developed  in  the  next  paragraph. 


Approximate  Solution  For  Shake  Test— Equation  (46)  tells  us:  (a)  when  the  system  is 
undamped,  i.e.,  when  [ii]  =  [0],  {a}and  {b}are  not  coupled  (the  solution  will  be  identical 
to  the  solution  of  Eq.  (10)  if  {a  }is  replaced  by  {q};  anu  (b)  when  the  system  is  damped,  or 
[i?  1  #  [0],{a}  and  {&}are  coupled.  Then  the  phase  factor  becomes  important,  and  a  set  of 
2 n  linear  simultaneous  equations  has  to  be  solved.  In  general,  the  semi-analytical  determina¬ 
tion  of  fixed-base  natural  frequencies  of  such  a  system  is  complicated. 


Since  numerous  small  viscous  damped  mechanical  structures  do  exist,  it  is  well  worth 
while  to  continue  the  investigation  of  these  systems.  In  order  to  establish  the  justification 
to  extend  the  developed  method  to  such  cases,  we  recall  Eq.  (46).  In  Eq.  (46),  one  sees  that 
when  damping  is  small,  the  real  part  of  the  response  {a  }  is  weakly  coupled  with  its  correspond¬ 
ing  imaginary  part  { b },  and  Eq.  (44)  approaches  Eq.  (10)  as  damping  decreases  towards  zero. 
Consequently,  the  solution  of  the  real  part  of  Eq.  (46),  i.e.,  {a},  also  approaches  the  responses 
{q  }  of  the  undamped  system  in  Eq.  (10),  which  is  real.  Therefore  for  a  slightly  damped  sys¬ 
tem,  the  solution  may  be  approximated  by  using  Eq.  (46)  for  shake-test  analysis.  To  be  more 
specific,  we  rewrite  Eq.  (46)  in  the  following  form 


[z]  j  -t 4*1  f{f}l 

u[R]\  [Z]  [{b}J  \{0}J  ’ 


(47) 


where  [Z]  is  tne  flexibility  matrix  of  the  very  same  system,  but  with  no  damping.  Then  the 
real  part  of  the  solution  of  Eq,  (47)  may  be  regarded  as  the  approximate  solution  of  a  slightly 
camped  system.  To  show  this,  we  suppose  that  the  system  under  study  is  not  excessively 
damped;  the  real  part  of  Eq.  (47)  may  be  approximated  as  the  solution  of  the  system.  In 
solving  for  the  responses  in  terms  of  applied  forces 
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(48) 


(49) 


{a}  =  {geos*} 

-  (r/1  +  u>2lZ)~HR ][Z]~l [R ]) _1  [Z H  {F }  (50) 

=  ([/I  +  U^fR])2)-1^]-^} 

where  the  matrix  u[Z]~l  is  the  mobility  matrix  HI.  When  [R 1  is  small,  the  following  ap¬ 
proximation  may  be  made: 

([/]  +  (l KHRl)2)-1  -  I/]-1  -  I/]  ,  (51) 

because  the  elements  of  the  matrix  [  ][R  ],  m,yry-ft ,  are  much  smaller  than  the  damping  ma¬ 

trix  elements  r,y ,  which  are  in  general  small.  Then  Eq.  (50)  becomes 

{geos*}  a  IZ]-1^}:.  (52) 


It  is  exactly  the  form  of  Eq.  (14),  except  that  the  responses  involve  a  phase  factor  cos  *. 
Therefore,  the  developed  method  can  be  applied  to  slightly  damped  system  as  well. 


EXAMPLE  PROBLEM  AND  DIGITAL  COMPUTER  SIMULATION 

A  specifically  designed  sample  is  in  order  to  illustrate  the  applicability  of  the  developed 
method.  Also  some  comparisons  were  made  to  bring  out  the  ambiguity  which  could  arise  due 
to  negligence.  A  simple,  but  still  general,  undamped  one-dimensional  system  of  nine  degrees 
of  freedom  as  shown  in  Fig.  1  will  be  treated.  The  equations  of  motion  for  free  vibration 
may  be  derived  from 

-  ~~  *  0  where  L  =  T-V  ,  (El) 

dt\dy,)  dy, 
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Fig.  1  *  Configuration  of  the  nine-degree-of-free- 
dom,  one-dimensional,  undamped  composite  struc¬ 
ture  for  distal  computer  simulation 


(yi  ->2)|z  +  y3  +  Lal^~  (y2~y3)l2 

]' 


*  +  |A#5[y5  +  ^J^(y4-h)]2*  I^e[y6 

+  \ M7yjr  +  |w8>72  +  |^9 ^8  +  L\~~  &  -y«>]2 


+  ^azia^-ye) 
1*2 


+  2‘Vl0 


-V9 


i>2  ~  ^10 
*2 


(ys~y9)j2  +  | 


Mny92  * 


(E2) 


v-  |Ki(yi-y4)2  +  |^2[y2  -  ^-^(y£+y4)j2 

+  |/f3|y3  -  L2LgL8  (ys  +y6)]2  +  |*4(y4-y7)2  +  l^ys-yg)2 
+  |^6(y6_y9)2  +  \Kw*  +  2 ^syg  +  \K*y$  •  (E3) 


For  this  particular  system,  the  kinetic  energy  T  is  coordinate  independent  and  the  potential 
energy  is  velocity  independent;  thus  the  Lagrangian  equation  has  the  form 


d/dT\  +  dV 
dt\dyij  +  dy; 


(E4) 
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By  differentiation  of  T  and  V,  the  equations  of  motion  are: 


Ml  +  j*1  +  +  *iyi '  *iy4  =  c’ 

•  Haf  •  -wfl*  • 

*  Kiy2  -  K2(^)yt  -  **(^)»  -  0, 

“3(£fcrs)(£:)*2  +  "3®)  ****  - s*[^rr^)n 

-  ^(^rrh '  °* 


)"-“.fi^©«--«."-'-(fifi)» 


+  ^2(^T~)  «■  *4  >4  +  *2(^11~)  y5  -  *4y7  -  0, 


-fifiHft)*.  *  [-(ft)'  • 

-  -fir*)"  -  ■■fir')"  •  ■•fir*)’" 

+  ^(^ZT3)  +  "  *5  y*+K3[^t^)  ye  *  =  o, 

„  (^-LsV'LsY,  A  L/M2  .  wlr.  „  (L2-Ls\ 


M6{~^r)[^)h  +  Me(&)  +  y«  -  Ka( 


^2  ~  ^8 


Ms  +  Mg 


'Ll  -  fV 


/Lj  -  Lg\  / Lg\  .. 


)  ys  +  K3[^lT*)  +  x6  y6  -  ^y9  =  o, 

—« 

2 

-)  y7  +  -  *4>4  +  (K4  +  K7)y7  =  o, 

h  +  M9(z?)  +Mio(~l7^)  y® 


y7  +  Mg 


_  K5y5  +  (Ke,  +  Ks)yg  _  0, 


(E5  Con’t.) 


18 


*10 


SI  AKDSKOf 


+  J#n 


>9  “  *6X6  +  (*6  +  Kgiyg  = 


0. 


(E5) 


We  tzy  a  solution  as  the  teal  part  of 


y,  =  ylVwf.  where  y,  =  foie*0-  (E6) 

Prom  straightforward  substitution  of  Eq.  (E6)  into  Eq.  (E5),  the  influence  matrix  \Z ]  can  be 
determined. 


M- 


The  elements  of  the  [Z 1  matrix  are  real,  in  this  case,  and  symmetric. 

zn  -  K,  -  •.-*(«,  *  . 

Z‘2  -  Z21  -  • 


Z13  ~  Z3l  ~  Z15  ~  Zbl  ~  Z16  ~  3si  “  Zn  ~  ^71  “  *18  "  *81  “  *19  “  *91  -  0  » 

*14  =  Z41  =  -  Kl  , 


Z  26 


Zn 
z  • 


*62  ~  *27  ~  *72  ~  *28  '  *82  ~  *29  ~  *92  ~  0  • 


*43  -  Z 37  -  Z73  -  Z>  -  Zh3  -  Z39  -  Z93  -  0  , 
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z“  -  z“  ‘  ■ 

zx-za-  . 


Z„  =  *,  .  ♦  K,  -  . 

z4S  -  zM  .  ^izUf  -  -*n(^Xz 0  • 
z*6  =  z6*  =  ^48  =  zM  =  z49  =  Z*  =  0  , 

■Z47  =  U74  =  -  K4  , 

Z“  ■  ^(^TT3)2  *  ^ 

z»-z«-^(^rs),-"I*«(^)(^). 


^57  “  ^75  “  ^59  "  ^95  '  0  , 

Z5S  =  Z8S  =  ~KS  , 

"ity* +  M\ 

Z67  ~  Z76  ~  Z68  ~  Z86  ~  0  , 

Z69  =  Z96  =  ~  K6  , 


z«  ’  *  K* 


Z77  -  Z4  ■»  JC7  ~  w2  jflfg  +  M,  p1  ^  j  1  , 

*»  =  z»’  -  ■ 

279  =  Z9"  =  0  , 

2  r 

Z88  =  tf5  +  K8  ~  W2  Af9^j  +  Af10|kio)' 
z69  -  z98  .  -  9>2.v,o(i^is)(g) , 


2 

Z99  =  K6  +  Kg  -  u2  zWio(|aj  +  3/ 


J 


(E7) 
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In  calculating  the  natural  frequencies,  the  numerical  values  of  the  physical  constants  are 


JT,  *  Jfs  ~  *3  *  1  . 


=  Jf5  =  J#S  =  2  . 
*7  -  1  . 

Jtfg  =  J#9  =  J#10  =  3  . 


*11  =  4  . 


Li  =  3  , 

Lj  =  L\  -  4  , 

Lz  =  I.9  =  2  , 

Ls  =  lsj  =  Lg  =  Lo  =  1  , 
=  2.5  , 

*1  =  K2  =  1  . 

K3  =  K4  =  1.5  , 

K5  =  K$  =  1  , 

K7  =  2  , 

K6  =  Kc  =  3  . 


The  natural  frequencies  of  the  total  system  are  calculated  by  setting 


The  results  are:; 


[z\:  =  o. 


u;  j  =  0.3631  , 
u-  2  =  0.4802  , 

uj-j  =  0. 1 801  . 

u.  [  =  0.9196  , 


4  =  1.0266 
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ul  =  1.2656  , 

=  1.4065  , 

=  1.4890  , 

=  2.1373  . 

Separate  cases  will  be  taken  into  consideration. 

Caae  1  —  General  Uastratioa  at  the  Applicability  of  the  fropoaed  Method  with  Forces 
Applied  at  the  Supporting  Points 

Define  the  subsystems: 

y\,  y2,  and  >3  span  the  subspace  of  the  equipment  subsystem  (Fig.  2), 

>4,  >5,  and  span  die  subspace  of  the  support  system, 
y7,  yg,  and  yg  span  the  subspace  of  the  base  subsystem. 


Fig.  2  -  Configuration  of  the  equipment  substructure 
for  illustration  cue  1,  with  forces  applied  at  its 
supporting  points  h  yg,  and  TG  and  at  paint  yz  on 
the  equipment 


k®— <§> 


K'i 


-® — 


♦*f5 


♦K 


In  carder  to  test  the  method,  the  fixed-base  natural  frequencies  of  the  equipment  subsystem 
are  first  calculated  by  setting: 


zu 

^12 

Z\z 

\[Ze)\  =  0  = 

z2\ 

Z22 

Z23 

£31 

ZZ2 

Zzz 

cjj  =  0.7816  , 
<4  =  0.8774, 
uj  =  1.2794  . 


The  results  are: 
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Then  apply  the  resonance  condition 

V4  ~  >5  =  ^6  =  0, 

with  forces  (3.  f4,  f$,  and  /g  applied,  one  at  a  time,  at  the  equipment  and  supporting  points 
y »,  y4.  y 5,  and  yg  respectively.  The  resonance  function  £(mj;-;u)  may  be  written  as: 


where 


J  mll 

m\% 

m15 

ml6 

jj 

1  m41 

m44 

m45 

m46 

f  m51 

mS4 

n,55 

m57  ! 

1 

i  m6i 

m64 

m65 

m66  j 

v(/ny;ui)  - 


m44 

m45 

"1 46 

mS4 

m55 

mS6 

i  m&i 

m65 

n*66 

_  _  %  _  toyi; 

"  '  “'ilzji  fj 


and  Z7  j  is  cofector  of  z/y  in  Z.  The  result  of  the  digital  computer  simulation  is  shown  in  Fig. 
3.  The  peaks  occur  right  at  the  frequencies  0.786, 0.8774,  and  1.2794. 


Case  2  —  Illustration  of  the  Actual  Physical  Location  of  the  Equipment  Subsystem  Being 
Immaterial 


In  this  case,  the  attempt  is  made  to  show  that  the  actual  physical  location  of  the  equip¬ 
ment  subsystem  is  immaterial,  wherever  response  devices  can  be  positively  secured,  and  meas¬ 
urement  can  be  taken.  The  systems  in  this  case  are  redefined  as:. 

v4>  y5,  and  y$  span  the  subspace  of  the  equipment  subsystem  (Fig.  4), 

yi,  v2,  V3,  y7,  ys,  and  yg  span  the  subspace  of  the  support  subsystem. 

The  rest  belong  to  the  base  subsystem. 
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Fig.  3  -  Resonance  function — frequency  piot  for  digital 
computer  simulation  of  iBustxatioa  caw  1 


Fig.  4  -  Configuration  of  the  equipment  substruc¬ 
ture  for  illustration  case  2,  with  forces  applied  at 
its  supporting  points  y  x ,  y2f  y$,  yz.yg,  andyg.and 
at  point  yg  on  the  equipment 


fueoucncy  m 


Again  solve  for  the  fixed-base  natural  frequencies  by: 


£44 

Z45 

Z46 

\[Ze]\  =  0  = 

Z54 

Z55 

Z56 

Z64 

Z65 

Z66 

The  following  listed  are  the  calculated  frequencies: 

=  1.0712  , 
u>2  =  1.3015  , 
uj  =  1.9682  . 
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The  resonance  conditions  are:. 

yi  =  yi  -  >3  =  >7  -  >8  =  y9  =  0  » 

with  forces  fu  f2.  /a.  h,  fs*  /9,and/6  applied,  eadi  at  a  time,  at  the  points  yi,  y2,  y3.  y7. 
y8,  y ; ,  and  y6  '^spectjveiy. 

The  resonance  function  is: 


Mwiy.-faj)  = 


"*66 

"*61 

"*62 

"*63 

"*67 

"*68 

"*16 

mn 

"*12 

"*13 

"*17 

"*18 

m  26 

"*21 

"*22 

"*23 

"*27 

"*  28 

"*36 

"*31 

"*32 

"*33 

"*37 

"*38 

"*76 

"*71 

"*72 

"*73 

"*77 

"*78 

"*86 

"*81 

"*82 

"*83 

"*87 

"*88 

"*96 

"*91 

"*92 

"*93 

"*97 

"*98 

"*11 

"*12 

"*13 

"*17 

"*18 

"*19 

"*21 

"*22 

"*23 

"*27 

"*28 

"*29 

"*31 

m  32 

"*33 

"*37 

"*38 

"*39 

"*71 

"*72 

"*73 

"*77 

"*78 

"*79 

"*81 

"*82 

"*83 

"*87 

"*88 

"*89 

"*91 

"*92 

"*93 

"*97 

"*98 

"*99 

The  mobility  elements  are  defined  as  before.  The  result  of  the  digital  computer  simulation 
shows  that  the  peaks  occurring  in  the  graphical  representation  truly  coincide  with  the  cal¬ 
culated  values,  which  are  shown  in  Fig.  5. 


Case  3  —  Illustrations  of  Frequency  Shift  Due  to  Change  of  Loading  Condition  and  of  Pos¬ 
sible  Ambiguity  Due  to  Coupled  (or  uncoupled)  Subsystems  within  the  Equipment 
Subsystem 

This  case  is  designed  with  the  intention  of  showing  that  any  detachment,  addition,  or 
rearrangement  of  loading  condition  within  the  equipment  subsystem  could  cause  a  large  shift 
in  the  fixed-base  natural  frequencies  of  that  subsystem.  (Actually  it  will  effect  all  the  natural 
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Fig.  5  -  Resonance  function— frequency  plot  for  digital 
computer  Emulation  of  illustration  case  2 


frequencies  of  the  total  system.  Here  emphasis  has  been  put  on  the  equipment  subsystem, 
because  it  is  the  main  concern  in  this  case.) 

Everything  remains  the  same  as  in  Case  1,  except  that  the  mass  Jlf2  which  couples  y* 
and  y2  is  removed  (Fig.  6).  The  fixed-base  natural  frequencies  calculated  by  setting  M2  =  0 
are: 


co*  =  1.000  , 
<o  2  =  0.786  , 
co  3  =  2.938  . 


When  Af2  is  set  equal  to  zero,  it  physically  decouples  yi  and  y2,  which  causes  the  existence  of 
uncoupled  subsystems  within  the  equipment  subsystem.  Caution  must  be  exercised  to  clear 
the  ambiguity  expected  by  the  theoretical  analysis.  In  order  to  reveal  all  the  fixed-base 
natural  frequencies  of  this  partially  decoupled  equipment  subsystem,  more  than  one  point  on 
the  equipment  must  be  shaken.  More  specifically,  we  start  with  the  influence  matrix  of  the 
equipment.  When  M2  =  0,  Z\2  =  Z2i  =  0, 


0 

0 

^22 

^23 

^23 

Z33 
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Fig.  6  -  Configuration  of  the  equipment  substructure 
for  illustration  case  3,  with  forces  applied  at  its 
supporting  points  y4 , ys ,  and  yg,  and  at  either  point 
y3  or  point  yj  on  the  equipment 


where 

Zn  =  *1  -  <o2 Mj  , 

^13  =  ^31  =  0  , 

Z22  =  *2  w2M3  (^ZT1)2  ' 

.at, 

2 

Z33  =  *3  -  «2  M +  M4 
When  the  shaking  point  on  the  equipment  is  at  y 3, 


co  2  0.786  , 

uc3  =  2.938  , 


ar.d  miss  the  one  u>\  =  1.000  (Fig.  7). 
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Fig.  7  -  Resonance  function- 
frequency  {riot  for  digital  com¬ 
puter  simulation  of  illustration 
case  3— no  force  applied  at  yj 


When  the  shaking  point  on  the  equipment  is  at  ?2,  the  result  will  be  the  same  as  shaking 
at  3%  because  in  the  resonance  function,  only  the  numerator  changes  from  Z22  to  Z33,  and 
it  will  not  affect  the  frequency  values  corresponding  to  the  peaks.  However,  if  the  shaking 
point  on  the  equipment  is  at  yj ,  then 


The  resonance  function  in  this  case  is: 

lll[Z‘]lll  IIZn(Z22Z33  -  z‘a)\\  IIZ11II 

The  corresponding  computer  simulation  is  shown  in  Fig.  8,  which  reveals  only  one  of  the 
fixed-base  natural  frequencies  coe  and  loses  the  other  two.  The  complete  information  has  to 
be  the  combination  of  these  two.  The  important  points  illustrated  here  are  the  following: 
first,  the  frequency  shift  due  to  change  of  loading  condition  (Fig.  9);  and  second,  the  possible 
ambiguity  which  may  arise  from  a  decoupled  equipment  subsystem.  To  the  second  point,  it 
appears  to  be  a  good  practice  to  take  one  or  two  more  shaking  points  on  the  equipment  for 
checking  purposes. 


(Z22  Z23 

zZ2  £33 

2 

w  \\z22ZZZ  ~  Z23^  _  “ 


Case  4  —  Illustration  of  the  Irrelevancy  of  the  Physical  Location  of  the  Shaking  Points  and 
its  Key  Rule 

This  case  is  designed  to  show  that  the  shaking  points  can  be  on  the  base  points.  With 
the  same  definition  of  the  subsystems  (or  subspaces)  in  Case  1,  the  influence  matrix  of  the 
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Fig.  8  -  Resonance  function— 
frequency  plot  for  digital  com¬ 
puter  simulation  of  intubation 
case  3— no  force  applied  at  yj 


Fig.  9  -  Resonance  function— frequency  plot  for  com¬ 
parison  of  frequency  shift— case  1  versus  case  3 


equipment  subsystem  remains  the  same.  The  theoretical  analysis  predicts  that  the  resonance 
function  is  independent  of  application  points  of  the  forces  although,  in  practice,  one  meas¬ 
ures  different  mobility  entities,  i.e.. 
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mu 

m41 

«51 

mgi 

^K-")  =  ITT 


"•17 

m18 

m19 

m47 

m*8 

m*9 

mS7 

mS8 

m59 

"*67 

m68 

m69 

m47 

m48 

m49 

m57 

mS9 

m67 

m68 

m69 

This  tinw>  the  shaking  points  are  at  yj,  yg,  and  yg .  The  computer  simulation  of  this  case  is 
shown  in  Fig.  10.  One  notices  that  they  exactly  coincide  with  the  result  obtained  by  apply¬ 
ing  forces  aty4,  y5,  and  y6-  The  same  result  will  be  achieved  by  applying  forces  at  the  com¬ 
bination  of  support  and  base  points.  Duplication  of  another  computer  simulation  will  not  be 
necessary.  The  key  rules  for  the  force  application  points  are: 

1.  Keep  the  shaking  points  away  from  the  equipment  subsystem. 

2.  The  number  of  shaking  points  is  equal  to  the  number  of  supporting  points 
plus  one. 


Fig.  0  -  Resonance  function— frequency  plot  for  digi¬ 
tal  computer  simulation  of  illustration  case  4— 
irrelevancy  of  the  physical  locations  of  the  shaking 
points 
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DISCUSSION 

The  theoretical  analysis  of  this  problem  is  rather  general,  as  long  as  the  dynamic  system 
under  consideration  is  undamped  or  slightly  damped  linear  time-invariant.  It  should  be  under¬ 
stood  that  a  reliable  result  in  actual  measurement  can  be  expected  only  if  the  test  is  being 
carried  out  according  to  this  condition.  The  proposed  mgrh*nir*l-rp«nnanrip  frequency- 
measurement  method  applied  to  a  multiple-support  equipment  system  is  quite  powerful, 
because  the  current  difficulties— for  instance,  the  requirement  of  farce  ratio  of  the  shakers 
(both  magnitude  and  phase),  physical  location  of  the  shakers,  and  size  of  the  object  to  be 
tested— are  removed,  or  their  limitations  are  lifted  partially  or  entirely.  The  method  suggests 
indirect  mobility  measurements  through  response  measurements  by  shake  test.  The  neces¬ 
sary  measuring  points  are  limited  to  only  a  few,  namely  the  support  points  and  one  or  two 
points  on  the  equipment  subsystem,  depending  on  each  individual  case.  Furthermore,  de¬ 
vices  and  techniques  to  perform  such  tests  are  currently  available,  and  results  can  be  acquired 
in  a  rather  routine  manner. 

The  theoretical  analysis  has  beer  proved  correct  by  computer  simulation.  A  vast  group 
of  mechanical  systems  and  structures  in  reality  can  be  treated  fay  this  analysis  as  if  they  are 
linear  time-invariant.  Therefore,  there  is  no  reason  why  it  should  not  be  feasible  to  determine 
the  fixed-base  natural  frequencies  semi-analytically.  We  have  been  talking  about  resonances, 
but  in  practice,  we  use  the  resonance  properties  only  to  deduce  the  fixed-base  natural  fre¬ 
quencies  without  actually  exciting  the  resonance  modes.  It  is  a  great  advantage  in  determin¬ 
ing  the  fixed-base  natural  frequencies,  because  in  reality  only  in  scarce  cases  may  one  excite 
a  resonance  mode  of  a  complex  subsystem  unless  “inner  resonance”  or  “beating”  exists.  It 
is  simple  to  illustrate:  in  the  analysis  of  the  text,  we  define  the  superstructure,  base  structure, 
and  substructure  on  an  arbitrary  basis.  There  is  no  difference  in  reality  between  the  super¬ 
structure  and  the  substructure,  both  physically  and  mathematically.  The  so-called  super¬ 
structure  requires  the  base  points  to  be  stationary  while  it  vibrates  at  its  fixed-base  natural 
frequency,  but  the  substructure  demands  otherwise;  therefore  they  interfere  with  each  other. 
That  is,  the  substructure  prevents  the  superstructure  from  vibrating  at  its  own  fixed-base 
natural  frequency.  There  is  only  one  condition  under  which  both  the  superstructure  and  the 
substructure  have  the  common  fixed-base  natural  frequency.  This  condition  is  what  is  called 
“inner  resonance,”  or  “beat”  phenomena.  It  is  also  said  that  the  system  is  degenerate. 


SUMMARY  AND  CONCLUSION 

It  is  realized  that  fixed-base  natural-frequency  analysis  of  a  mechanical  system  is  of  im¬ 
portance  in  the  dynamic  design  and  dynamic  behavior  study.  A  resonance  function  has  been 
derived  for  an  n  degree-of-freedom  system  to  measure  the  fixed-base  natural  frequencies  of 
the  equipment  subsystem.  The  analysis  is  based  on  a  linear  time-invariant  system. 

The  resonance  function  is  an  explicit  function  of  mobility  entities,  which  is  composed 
of  intrinsic  properties  of  the  total  system  with  vibrational  frequency  as  a  parameter.  It  is 
clear  that  as  long  as  the  dynamic  distortion,  or  disturbance,  remains  within  the  linear  limit, 
this  function  remains  invariant  with  respect  to  its  input,  output,  and  time  as  well.  This  in¬ 
variant  feature,  besides  the  advantageous  facts  pointed  out  in  the  discussion,  makes  the  pro¬ 
posed  method  more  desirable  than  currently  available  methods. 
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